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SINGULARITIES OF THE RADON TRANSFORM 



A. G. Ramm and A. I. Zaslavsky 

Abstract. Singularities of the Radon transform of a piecewise smooth function 
f(x), x S R n , n > 2, are calculated. If the singularities of the Radon transform 
are known, then the equations of the surfaces of discontinuity of f(x) are calculated 
by applying the Legendre transform to the functions, which appear in the equations 
of the discontinuity surfaces of the Radon transform of f(x); examples are given. 
Numerical aspects of the problem of finding discontinuities of f(x), given the discon- 
tinuities of its Radon transform, are discussed. 



I. Introduction 

Let f{x) be a compactly supported function, D be its support, and T — dD be 
a union of finitely many C°° hypersurfaces Ti, . . . ,T S in general position, each of 
which can be written in local coordinates as 

x n =g(x'), x' = (xi,...,x n -i),n > 2, 

where g[x') G C°° , f(x) € C°°(D), f(x)\ r > c > 0. The discontinuity surface of 
f(x) is F, the boundary of D. We assume that the rank of the Hessian gij(x) := 
d 2 g l 'dxidxj is constant on each of Tj, 1 < j < s. 

Define the Radon transform (RT) of f(x) by the usual formula [GGV] f(p, a) = 
L„ f(x)5(p— a-x) dx, where 5 is the delta-function. It is well known that f(Xp, Xa) = 
|A| _1 f(p 7 a), A e R 1 , A 7^ 0. Consider the integral 



(1) R(p,a;f) := / f(x)p(dx), 

Ji ap 

where l ap is the plane a-x — p = 0, a £ R", p £ R 1 , and pi(dx) is the Lebesgue mea- 
sure on l ap . One has R(p,a;f) — f(p/\a\, a ), a :— so that R(p,a;f) — 
\a\f(p,a), M = K + --- + ^) 1/2 - 

The problems we are interested in are: (PI) Find the singularities of R(p, a; /); 
and (P2) Find the surface T of discontinuity of f(x) given the singularities of 
R(p,a;f). 

No results concerning (P2) were known. In [N] one can find an estimate of the 
norm of f(p, a) in Sobolev spaces. This result does not give information about (PI) 
and (P2). In [P] there is a result given without proof, which has a relation to (PI). 
Our result is more general. In [Q] it is mentioned that the values (a: p), such that 
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l ap is tangent to T, play a special role. This observation is made quantitative in 
our Theorem 1. Our results are useful for inversion of incomplete tomographic data 
[R2]. 

The basic results are formulated in §11. They give solutions of the problems 
(PI) and (P2). Actually, more general problems are solved; particularly, finite 
smoothness of f(x) and T is allowed, the role of the intersections of Tj in the study 
of the singularities of R(p, a; /) is clarified, etc. In §111 proofs are sketched. In §IV 
examples are given. In §V numerical aspects of problem (P2) are discussed. 

We conclude this introduction by an outline of our ideas. First, we describe the 
behavior of R(p, a; f) in a neighborhood of the set Q / which is the set of singularities 
of R(p, a; /). Second, we prove that, in general, there is an equation of the set Qf 
which is of the form q = h{fi), [3 e i? n_1 , so that Qf is a hypersurface. Third, 
we prove that the function g(x') (in the equation of T) is the Legendre transform 
of the function h(/3) (in the equation of Qf). Fourth, we describe some geometric 
properties of Qf. 

Our results give a theoretical basis for the solution of the practically important 
problem in nondestructive evaluation and remote sensing, the problem of finding 
the discontinuities of a function from the knowledge of its RT. 

II. Formulation of the results 

The RT, defined by formula (1), is a function on the projective space MF n , and 
we take i?(l,0;/) := for compactly supported /. Let Qf denote the set of the 
points (a : p) in this projective space, which correspond to the planes l ap tangent 
to T = dD. We say that l ap is tangent to T at a point x € B m := HJli Fj, if l ap is 
not transversal to B m at the point x. 

1. Our first result is the following theorem in which the description of the singu- 
larities of R(p, a; f) is given. 

Let l ap be tangent to T at the point x. We claim that if a is generic, then the 
set Qf is a smooth hypersurface in a neighborhood U of (a : p) . If A is a symmetric 
matrix with real-valued entries, then its inertia index (inerdex) is defined to be the 
number of its negative eigenvalues. Consider first the case when V consists of one 
surface. 

Theorem 1. There exists an equation ((a: p) = 0, V( ^ in U , which defines 
Q f in U, and two C°° functions r\ and r 2 in U such that 



(2) R(p,a;f) 



>(n-l)/2 . .j. . 

ri +r2, ij In is even, 

C ( "- 1)/2 (ln|C|)ri +r 2 , if In is odd. 



Here I is the inerdex of the matrix Zkj, where Zkj is the Hessian of the function 
z = (a ■ x — p)/\a\ on T at the point x and z + = max(z, 0). 

If x € B m and (a: p) is generic, then the following result holds. 

Theorem 1'. There exists ((a: p), VC ^ in U , such that the equation ((a: p) = 
is the equation of Qf in U, and two C°° functions n and r 2 in U, such that 



(2') R(p,a;f) 



/ -(n+m—2)/2 . » T , i\ • 

C\_ ri+r2 1 if 1 [n + m — 1) is even, 

C ( " +m - 2)/2 (ln|C|)n +r 2 , if /(n + m-1) is odd. 
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In [RZ1] the constant r\(a: p) is calculated. In [RZ2] this result is used for 
a derivation of the asymptotics of the Fourier transform of a piecewise smooth 
function. 

2. Let us define the Legendre transform of a function g(y), y E in a neighbor- 
hood Uy of a point y at which the matrix gij(y) ■= d 2 g / dyidyj is nondegenerate, 
i.e., detgij(y) =/= in Uy. Define Lg := h{[3) := (3-y — g(y), where the dot stands for 
the inner product and y = y(j3) is the unique solution of the equation (3 = Vg(y) 
in a neighborhood U-p- of the point (3 — Wg(y). One can prove that if g E C l (Uy), 
I > 2, and dct g^y) ^ in U v , then h(j3) E C l (Up). 

It is known that under our assumptions Lh — g{y), i.e., the Legendre transform is 
involutive: g{y) = (3-y~h{(3), where (3 = (3{y) is the unique solution to the equation 
y = Vh(P), (3 E U-0. One can prove that det hij ((3) ^ in U-p if det gij(y) ^ in Uy', 
moreover, the matrix hij {(3) is inverse to gij(y), where (3 = (3{y). Recall that F is 
a union of hypersurfaces Tj, 1 < j < s, T\, . . . , T s are C°° and in general position. 
Denote B m := ri v .. ;TO the set of (a: p) € MP„ such that l ap is tangent to B m . The 
set B m C KP„ may not be a hypersurface (see Theorem 3); however, as Theorem 
1' claims, it is indeed a smooth hypersurface outside a set of (n — l)-dimensional 
Lebesgue's measure zero. 

3. Our second result gives the relation between the discontinuity surfaces for R(p, a; f) 
and those for f(x); namely, the function g(x') in the local equation o/F, x n = g(x'), 
is the Legendre transform of the function h{(3) which gives the equation of Qf, 

q = h(J3). 

Assume that q = h([3), (3 € U-p, where Up is a neighborhood of a point (3, 
q = h(j3), and det/i y (/3) ^ in Up, where hij := d 2 h/ dfrdfij. Let x' = Vh(f3). 

Theorem 2. // h(J3) E C l {Up), I > 2, then Lh = g(x'), and g(x') E C l (U w >). 

This result allows one to recover the surfaces of discontinuity of f(x) given the 
surfaces of discontinuity of R(p, a; f). 

4. Examples show that the Legendre transform h((3) = Lg of a function g(x'), 
x' E K n_1 , may have domain of definition of dimension less than n— 1. Since Qf is 
a union of several varieties of codimension one in MP„ (called components below) , 
the question arises: which of the components of Qf and which of their intersections 
provide, after applying the generalized Legendre transform defined in [RZ1], parts 
of r = dD which have codimension one in K™. The answer is given in Theorem 
3. This theorem describes Qf in terms of differential geometry of T. Recall that 
the principal curvatures of a hypersurface S C R™, which is the graph of a function 
x n = g(x'), are the eigenvalues of the matrix (g lj )- (% ■ (1 + gf)~ 1/2 , 
gi = dg/dxi. One can prove that if k, k> 1, principal curvatures of a hypersurface 
S vanish identically, then for every point P E S there exists an affine k- dimensional 
space Lp such that P E Lp C S . 

Theorem 3. (a) Assume that B m is nonempty. Then m principal curvatures of 
B m vanish identically; 

(b) // k principal curvatures of Ti vanish identically, then T\ has codimension 
k + 1 in MP„. 

Every point ofTi is a vertex of a cone K, which belongs to T\j, where T x f\Tj ^ 
0. The directrix of K is (k~ 1) -dimensional, and this directrix can be described as 
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follows: Take an arbitrary point P £ Ti, and let Lk(P) C Ti be a k- dimensional 
affine space containing P , which exists since k principal curvatures of T\ vanish 
identically. Let dp :— {{a: p): l ap be tangent to Tj at the points of Lk{P) H Tj}, 
and let l aoPo be tangent to Ti at the point P. The vertex of K is the point (a : p ). 
The directrix of K is the set dp. 

The set Qf is a union of the sets T il ... ik , Qf — {jT il ... ik where the union is taken 
over all combinations of indices 1 < ik < s. Theorem 3 gives a recipe to select the 
components of Qf which yield after the Legendre transform the components of T 
of codimension 1, i.e., hypersurfaces Tj which are parts of F, T = \Jj =1 Tj. Note 
that if a component of Qf has some principal curvatures vanishing identically, then 
its preimage in R™ has codimension greater than one. Therefore, if one wishes 
to recover hypersurface-type components of T, then one should apply the Legendre 
transform to those components of Qf, which do not have principal curvatures which 
vanish identically. Those hypersurfaces Tj which have identically vanishing princi- 
pal curvatures are reconstructed by applying the generalized Legendre transform, 
which was introduced in [RZ1], to high-codimension parts of Qf described in The- 
orem 3(b). The generalized Legendre transform was applied in [Z] to the study of 
dual varieties in algebraic geometry. 

It is well known that the Radon transform may be considered as a Fourier integral 
operator, so it makes sense to study its action on the wave front set of /. In [RZ1] 
we study a relation of the wave front of / and the set Qf. 



We sketch the proofs in the simplest case m = 1, n = 2, but the ideas are similar 
in the general case. 

First we prove that if D C C°° and / G C°°, then R{p, a; f) e C°° on the set 
Vf := WP n \Qf. Thus, the singularities of / are in the set Qf. Second, we prove 
that, generically, Qf is a C°° hypersurface in RP„ and find the equation of this 
hypersurface. 

Third, we prove that there exists a neighborhood U of a generic point (a : p) and 
an equation ((a: p) = 0, V( ^ in U, such that (2) holds. 

(a) Let us start with the second claim and prove also Theorem 2 for n > 2. Let 
Q'S-p=0bca tangent plane l ap to T at a point icT. Assume that a n ^ 0, 
and write x n = [3 ■ x' - q, fa := -ai/a n , q := -p/on, x' = (xi, . . . ,x n -\). Let 
x n = g(x') be the equation of T in a neighborhood U of x, and det gij(x') ^ 0. 
Then Vg{x') =/3,q = [3-x'-g{x'). Thus q = h(/3) := Lg. The equation q = h(}3) is 
the equation of Qf in the inhomogeneous coordinates (j3, q) . One can prove that if 
q G C S (U), s > 2, and det gij{x') ^ 0, then h e C S (U), where U is a neighborhood 
of the point (f3,q), Vg(x') = [3, q — [3 ■ x' — g(x'). Since L is involutive, g = Lh. 
Theorem 2 is proved. 

(b) Let us prove the first claim for n = 2. Assume that (a: p) £ Vf, i.e., l ap is 
not tangential to T. Write R(p, a; f) as 



where ttj := a^q, (3) are the points of intersection of l ap with T. The integral J is a 
sum of the integrals over the intervals (oi, b), (b, c), (c, a 2 ), where a\ < b < c < a 2 



III. Proofs of Theorems 1 and 2 
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and b, c do not depend on q, (3. Obviously the integral over (6, c) is a C l function of 
(3 and q if / e C l , I > 0. The integrals over (ai, b) and (c, 02) are treated similarly. 

Let us prove that the integral over (ai,b) is C l function of q,j3 if T,f e C l , 
I > 2, and l ap is transversal to T, that is, (3 ^ g'(ai)- It is sufficient to prove that 
ai(g,/3) <E C l . The function a\(q,(3) is the root of the equation q = /3ai — g{a\). 
By the transversality condition (3 — g'{ai) 7^ 0. Thus, the implicit function theorem 
implies that the root ai(q, (3) G C l if g <G C l . The first claim is proved. 

(c) Let us prove the last claim. Let (57: p) e Qf and ((3,q) be the corresponding 
nonhomogeneous coordinates. For a generic (a: p) the condition g"{xi) ^ follows 
from the equation g'(xi) — [3 and Sard's theorem. We can assume therefore that 
g"ix\) 7^ 0. Consequently, the point x\ is a Morse-type (nondegenerate) critical 
point of the function z := a- x — p on Fn U , i.e., of the function —(3x\ + g{x\) + q. 
The part of integral (1) taken over the complement to U is a C°°-function of 
(a: p) according to (6). It gives in formula (2). By the Morse lemma, there 
are coordinates ui,U2 such that the equation of T in these coordinates is u\ = 0, 
the region D n U is described by the inequality u\ > 0, and z = u\ + u\ in U. To 
study the singularity of R{p, a] /), take a curve 7 which intersects Q/ transversally, 
for instance, 7 = {(a: p): a — 57}. Parameter p gives the position of a point on 7. 
On lap one has 5 • x — p = and z = a ■ x — p, so z = p — p on lap- Thus, z can 
be used as a parameter which determines the position of a point on 7; therefore, 
the domain of integration in (1) can be described by the inequality u\ > and the 
equation z — ui — u| = 0. Thus, z — vi\ = u\, so — z]/ 2 < u 2 < z]/ 2 since z = z + > 
in the integration region. We have 

1/2 

R(p,a;.f) = / f{x)n{dx)= \ /i(iti, u 2 )m{du) = / f 2 (u 2 , z + ) du 2 , 

where f 2 {u 2l z) is a C°°-function, I is the curve given by the equation z — ui — u 2 = Q 
and iti > 0, n\(du) comes from /i(dx) via the Morse lemma change of variables, 
and the last integral comes after an elimination of u\. From this formula one 
derives (2). Indeed, write f 2 {u 2 , z) as a sum of even / e and odd / functions of u 2 , 
fe(u 2 ,z) e C"°°, f {u 2 ,z) e C"°°. Then the integral 

1/2 1/2 

r+ 1/2 r+ 

/ f c {u 2 , z+) du 2 = z+ n and / fo(u 2 , z+) du 2 = 0, 

where ri e C°°. The function r2 in formula (2) vanishes if T is strictly convex so 
that lap intersects T at two points only. 

IV. Examples 

1. Let f(x) = 1, \x\ < a, f(x) = 0, \x\ > a, x £ R n , n > 2, f(p, a ) = -p 2 , 
a° = a|o;| _1 . Thus p 2 /\a\ 2 = a 2 is the equation of Qf. In ((3, q) coordinates the 
equation of Q f is q = ±ay/l + /3 2 , (3 G K" -1 . Thus h{(3) = ±a^l + (3 2 . By 
Theorem 2 the equation x n — g(x') of the surface of discontinuity of f(x) is given 
by g(x') = Lh = =fV» 2 — x/ 2 . The equation x n = ±\ / a 2 — x' 2 defines the sphere 
\x\ = a. 

2. Let f(x) = 1, b < \x\ < a, f(x) = 0, \x\ < b or |x| > a, < b < a, n > 2. Then 
f(p, a ) = 2^a 2 -p 2 , b<p<a- f(p, a ) = 2(^ 2 ^ 2 - ^W^p 2 ), < p < 6; 
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f(p,a°) — 0, p > a, \a°\ = 1. Thus p 2 = \a\ 2 a 2 and p 2 = \a\ 2 b 2 are the equations 
of Qf. Taking Legendre's transform yields the surfaces \x\ — a and \x\ = b of 
discontinuity of f{x). 

3. Consider f{x) = outside of the region D bounded by T, where T is the union 
of the curves x 2 = and x 2 = x\ — 1, and let f(x) = 1, x £ D. The R(p, a; /) 
is a function whose support is bounded by the curves q = (3, q — —(3 from below, 
Q = \fl 2 + 1 in the interval — 2 < (3 < 2, and q = j3, q = —(3 for \j3\ > 2 from above. 
One can check that on the lines q = ±[3, — oo < [3 < oo, the function R(p, a; f) has 
a singularity of the type \z\ and on the parabola q = \(3 2 + 1 it has the singularity 

of the type z]/ 2 . Applying Legendre's transform first to the function q = j(3 2 + 1, 
—2 < (3 < 2, yields the parabola x 2 — x\ — 1, — 1 < X\ < 1; and secondly, applying 
it to the functions q = ±(3 yields two points x x — ±1, x 2 = 0. By Theorem 3, the 
straignt line joining these two points also belongs to T. Thus T is recovered. 

V. Numerical aspects 

The RT of f(x) is usually given with an error. Hence, the first numerical problem 
is to calculate the function h{(3) which gives the equation of the set Qf of the 
singularities of RT given the noisy measurements of the RT. The second numerical 
problem is to calculate Lh = g(x'). Calculation of the Legendre transform of a 
function h((3) known with errors is a well-posed problem, at least in the case when 
det gij(x') 7^ 0. It is proved in [RSZ] that if a function gs(x') given such that 
\gs(x') — g{x')\ < S, gs(x') is not necessarily in C 2 but is continuous, then one can 
calculate Lg with the accuracy O(S) as 5 — > 0. This means that a stable method is 
given in [RSZ] for calculating the Legendre transform of noisy data. See also [R5]. 
Our result in part 3 of §11 has an interesting connection with the envelopes theory 
[T, Zl]. 
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